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ABSTRACT 


A theoretical prediction to determine the frequency response of 
a:cascaded remote pressure transducer system has been developed. 

This prediction was used to find a theoretical configuration which would 
produce a fairly flat frequency response (1.00 + 0.05) overa 100 hertz 
range. Such a configuration would allow cross-correlation of an unsteady, 
random, low-speed flow field by simple analog methods. 

The theoretical prediction did not hold true for configurations of 
three serially connected tubes, but the parameters, namely, transducer 
volume, tube length, and tube diameter, proved the same trends experi- 
mentally that were predicted theoretically. By observing these trends, 


a> solution of 1.00 £0.09 overa 95 hertz range was obtained. 
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Because of the increasing interest in measuring oscillatory and 
random unsteady flow fields on lightweight aeroelastic models, a need 
has gradually developed over the past several years for a direct measuring 
remote pressure transducer system to allow rapid and economical deter- 
mination of auto- and cross-correlation type information plus the attend- 
ant spectral distribution. Such a system has been developed by Bergh 
(2 at the National Aeronautical and Astronautical Research Institute 
(N.L.R.) in Amsterdam, Holland, utilizing precise lengths of tubing to 
connect the remote pressure transducers with static orifices located on 
the model surfaces. Johnson [2] has confirmed Bergh's one tube theory 
and Allen (3) utilized a two tube configuration: both of these systems 


required compensation for the system's transfer function to be included 


z= 
‘ 


in the data reduction process. 

Should a configuration be found which would produce a flat fre- 
quency response overa specified range, cross-correlation of an unsteady, 
random signal could be achieved directly through the use of an analog 
circuit. Since a direct analog process does not readily allow for cor- 
rection of the system's transfer function, the primary purpose of the 
subject program was to stg a useable flat frequency response of 
1.0 £0.05 over a 0 - 100 hertz frequency range. This would allow 
direct evaluation of the unsteady aerodynamics in many low-speed flow 


Prellas ot interest. 


1] 








II. EXPERIMENTAL PROCEDURE 


A. GENERAL 
All experimentation was performed either in the Aero-Mechanics 
Laboratory, Room 022, Halligan Hall, or with the IBM 360-67 digital 


computer available at the U. S. Naval Postgraduate School. 


B. THEORETICAL BACKGROUND 

Bergh [3] has developed a theoretical prediction for the pressure 
ratio of not only a single tube remote pressure transducer system 
(Figure la) but also for a cascaded tube system employing N tubes 
such as shown in Figure lb. The solution is based on several assump- 
tions; the perturbations are of a sinusoidal nature imposed on the average 
values of the parameters, fluid flow in the tubes is assumed laminar, 
and the tube lengths are much larger than their aor By introducing 
these assumptions, the governing equations, i.e., the Navier-Stokes, 
continuity, energy, and state equations are reduced to a set of partial 
differential equations which may be solved subject to satisfying the 
boundary conditions. The mass flow rates of the air in the system 
between any two adjacent tubes must be equivalent. This relation pro- 


vides an equation which reduces to the following form: 


Ps wry W; ls Tg@™;) 

oe. Cc -] - 5 ae eee ee Minato WZ. o fe 
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where: 


/ ( 
p.= 2 yh aie FL Teles Pe) 7% 
3 Gee T2 («3 ) x a («5 PL Mz) 





a aad Ae Rs ( Ss) 


When a single tube configuration (see Figure la) is used, the entire 
third term of equation (1) disappears, simplifying the equation to the 


form below: 


P53 Opee ye (Vig Lene doe.) =| 
—— =]cosn ($3b5)+ 4 sin (6, L; 
Pye Ma)t aay Vts O35 Tix;) (#3 ) 








The subscript j refers to the tubes being analyzed in the N tube system. 
From equation (1) it is seen that the pressure ratio prediction is a func- 
tion of tubes j-1 and j+l as well as tube j, which in turn presents diffi- 
culties when optimizing the problem for the desired result. The solution 
complexity is not obvious at first; however, a closer inspection bears 
out that all the Bessel functions are functions of the imaginary parameter, 
oe which are sometimes called Kelvin functions (4) and hence, the 
solution is a complex quantity. The complex pressure ratio specifies 

not only the pressure ratio magnitude but also the phase shift. The 
equation also contains two uncertain quantities; the increase in trans- 
ducer volume due to diaphragm deflection, and the polytropic gas 
constant, k. Bergh proved both of these terms to have only minor effects 


on the solutions to equations (1) and (2) and hence adiabatic flow was 


assumed while transducer volume changes were ignored. The system 
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can be represented by the acoustic analogy provided by Olson [5] for 

a multi-degree of freedom system. The analogy develops equations of 
motion for a sound wave of pressure p travelling in a tube and impinging 
oma cavity volume, with the resistance represented by the viscosity and 
the capacitance term dependent on the cavity volume. Similarly, Bergh's 
solution predicts the complex pressure ratio with the resistance and 
capacitance again represented by fluid viscosity and cavity volume, 


respectively. 


oF COMPUTER TECHNIQUE 

‘With the aid of Bergh's equations, the pressure ratio at a certain 
frequency can now be predicted. However, because the pressure ratio 
of the next higher tube is part of the predicted ratio and the tube size 
and cavity volume interact in an irregular manner, the total solution 
does not lend itself for routine optimizing techniques to produce a flat 
frequency response. Therefore the equations were programmed for a 
digital computer to facilitate pressure response prediction using pro- 
gram-operator chosen configurations. The program was designed for 
use on the IBM timesharing system but can also be used for batch pro- 
cessing, and contained a program option allowing the selection of 
either a one-, two-, or three-tube configuration. Each tube stage ina 
multiple tube configuration is specified by three parameters: transducer 


volume, tube length, and tube diameter. The program was verified for 


proper operation by comparing it with both a hand-worked calculation 


US 








for a one-tube configuration and Bergh's [1) results for both one- and 


two-tube configurations. 


1am EQUIPMENT 
i Pressure Transducers 

The two pressure transducers were designed by 
Professor L. V. Schmidt [6) , and consisted of a Bentley Detector 
System mounted behind a 0.003-inch thick annealed brass diaphragm 
which was inserted in an aluminum housing. The diaphragm thickness 
was selected to give optimum sensitivity ina H 0.75 psi range. The 
detector system consisted of a Bentley 316 AJW proximity detector 
inserted into a brass diaphragm holder. The proximity detector was a 
reluctance gage which sensed change in the air gap between the reluct- 
ance coil and a metal target, in this case the brass diaphragm. The 
diaphragm deflection was directly proportional to an applied pressure 
perturbation, and the proximity detector output was linear with respect 
to the diaphragm deflection; therefore, the output was a relative measure- 
ment of the pressure applied to the diaphragm. 

Le Calibration Chamber 

The calibration chamber was the same as used by Johnson 
[2] and consisted of an aluminum block with a 1.5-inch diameter 
cylindrical chamber 6 inches long. The chamber ends were designed to 
accept either the transducer holders or the adapter plates which connected 
the chamber to a remote transducer through the transmitting tubes. The 


Chamber had a pressure port to allow the mean static pressure to be 


16 








varied from an external source, and also a threaded orifice designed to 
accept either a solid plug or a University Sound ID-75 acoustical driver 
unit. Both orifices were located at the chamber's midpoint so that 
unsteady pulses were introduced symmetrically to both transducer 
systems, namely the reference transducer at one end and the complete 


transducer system at the other end. 


le CALIBRATION 
ie General 
All calibration procedures were identical to those performed 
by Johnson ]2) and Allen {3] 
Ze Proximeter Calibration 
The Bentley 316 AJW proximity detectors were usually cali- 
brated in order to determine the voltage output with respect to the air 
gap between the reluctance coil face and a metal target. By mounting 
a metal target on a micrometer stem, the target could be backed off 
known distances from the reluctance coil and the output monitored on 
a D.C. voltmeter. The Bentley D252 carrier amplifiers have a variable 
amplification capability which may be adjusted until the proximity 
detector output varies proportionately with the target distance. Because 
of the large apparent sensitivity of the proximity detectors, this pro- 
cedure was foregone and the detectors were calibrated with respect to 


the pressure differential across the diaphragm during the static calibra- 


tion procedures. The proximity detectors were inserted into the brass 


a7 








diaphragm holders until their output was -3.0 volts D.C.: this was the 
range where transducer output and diaphragm deflection varied linearly. 
4 Static Calibration 

The transducer output was connected to a Bentley D252 
carrier amplifier powered by a -18.0 volt D.C. regulated power supply. 
By inputing the carrier amplifier output simultaneously with a biasing 
voltage into an Astro-Data D885 differential amplifier, the differential 
amplifier output, monitored ona D.C. digital voltmeter, was adjusted 
to read zero for zero pressure differential across the diaghaen. A 
U-tube manometer and hand crank pump conveniently allowed known 
pressure differentials to be introduced into the calibration chamber for 
static sensitivity determination. The complete instrumentation diagram 
is shown in Figure 2. Each proximity detector was calibrated during 
this stage by introducing positive and negative pressure differentials 
of equal magnitudes into the calibration chamber and adjusting each 
carrier amplifier until the output was approximately the same for equal | 
but opposite pressure differentials over a -0.8 to +0.8 psi range. The 
static calibration curve is shown by Figure 3 and bears out the linear 
sensitivity variation as measured at the carrier amplifier output of 1.37 
and 1.04 volts/psi for transducers 1 and 2, respectively, over a 40.8 
psig pressure range. 

4, Dynamic Calibration 
For the dynamic calibration, a sinusoidal perturbation was 


introduced into the chamber by closing the static pressure port and 
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inserting an ID 75 University Sound acoustical driver unit driven by a 
Wavetek 114 signal generator which was amplified by a power amplifier. 
The acoustic driver unit was adjusted to operate at an input level of 
13.5V RMS which corresponds to an oscillating pressure amplitude in 
the calibration chamber cavity of approximately 0.13 psi at 120 hertz. 
The carrier amplifier output was monitored with both a true RMS meter 
and a precision phasemeter, and visually monitored on a dual beam 
oscilloscope. The dynamic calibration instrumentation diagram is shown 
in Figure 4. The frequency was varied from 10 - 120 hertz in S-hertz 
increments. Since the disturbance was introduced symmetrically into 
the chamber, there was no appreciable phase shift in oscillating 
measures at either end of the chamber; however, the sensitivity of 


transducer #2 was 0.787 of transducer #1 at every frequency increment. 


Bs EXPERIMENTATION 

The data runs were made with the same system as was used for 
the dynamic calibration procedure with one exception; transducer #2 
was remotely located from the calibration chamber through the use of 
nylon and stainless steel tubing. The brass diaphragm had a centered 
circular wrinkle which allowed brass shims to be inserted between the 
diaphragm and the plexiglass spacer in the transducer housing without 
hindering movement of the diaphragm. The wrinkle volume was estimated 
as 0.00211 in.° and the shims had an inside diameter of 0.523 in.: 
this allowed any desired transducer volume to be constructed. Frequen- 


cies of 10-120 hertz were investigated in 5-hertz increments for 
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different one-, two-, and three-tube configurations . The ratio of trans- 
ducer output #2 divided by transducer output #1 was calculated and, 
because the sensitivity of transducer #2 was determined during dynamic 
calibration of transducer #1 as 0.787, the quotient was divided by 

0.787 to normalize the output of transducer #2 with respect to transducer 
#1. This determined the frequency response (transfer function) of a given 


configuration. 
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II. DISCUSSION AND RESULTS 


The computerized solution of Bergh's [1] prediction allowed 
quick and easy determination of the frequency response for a one-, two- 
or three-tube system, and could be easily modified for a greater number 
of tubes. 

The one-tube solution generally followed the typical organ pipe 
solutions which produced sharp resonance peaks at many frequencies. 
When the resonance peaks did not occur, a heavily damped solution 
appeared. A solution was sought that would reduce the resonance peaks 
or raise the damped solution to the desired goal (1.00 + 0.05); however, 
no such solution could be found fora 0-100 hertz range. The data runs 
did confirm Bergh's theoretical predictions for both pressure ratio and 
phase shift for a one-tube configuration, as can be seen by Figure 5. 
Therefore, a two-tube ee was examined in hopes that the pressure 
of an additional tube would counter the undesirable characteristics of 
the first tube. —_ a configuration was found with the desired frequency 
response of 1.00 + 0.05 from 10 to 100 hertz as shown by Figure 6. 
However, it was felt that the first tube in the system had a diameter 
that would result in the static orifice on the body being large enough to 
perturb the flow in the region of static pressure measurement. The two- 
tube configuration also verified Bergh's theoretical prediction of phase 
shift, but the pressure ratio was consistently lower (6%) than predicted 


for frequencies greater than 40 hertz (see Figure 7). A three-tube 


24 





PRESoURESRATIO, Po/P, 


PHASE SHIFT, DEGREES 


_.... THEORETICAL |: © 
© EXPERIMENTAL 


—.= ITHEORET ICAL 
© EXPERIMENTAL 





TRANSDUCER VOLUME (CU. IN.) OS LSS 


UIE WIRING Oe 8 IN) 23) 10 
TUBE DIAMETER (IN.) Or 07'S 
PiGU RE so 


SINGLE TUBE CONFIGURATION 


ZS 





SNOILNIOS IVOILIVOIHL AIYISIAG 
9 FUNnoId 
ZLYSH ‘AONSTNDOAYSA 
Oeil OOT Og 09 Ov OZ O 
| | I 
Vera LjHO< 0) > serrA0) 8 (0, S6G6.0 110°0 (°NI) “VIG 3ENL 


oe? G°6¢ €°O Oc Ce Cue evel) eldest 
GEZO'O 0°O 0°O GoGo 0 © © (- NI) ° TOA *SNVYUL 

















Ga Gs ee Ca ieee: 
Jan L Sqn 1 Sdn lL JEN L Jan L 
eR pepe el We) a ee 





NOILVYNOIANOD AGNL © ———-— 
NOTIVGnO Tange sad 





ae 


68) 
© 


O) 
© 


@ 
oH 


— 
e 
et 


Ta/Vq 9 la/€dq ‘oILva aYunssaud 


26 





PRESSURE RATIO, P3/P, 


PHASE SHIFT, DEGREES 










SS Weel Geb 
© EXPERIMENTAL 


0 20 40 60 80 100 120 
FREQUENCY, HERTZ 


O ae 40 60 80 100 N26 


—— THEORETICAL 
© EXPERIMENTAL 


100 





TRANS. VOL. (IN.?) 0.0 0.0238 

TUBE LENGTH (IN.) me (6 Dey 

TUBE DIA. (IN.) 0.078 0.0465 
FIGURE 7 


TWO TUBE CONFIGURATION 


27 





configuration was then sought in order to reduce the diameter of the first 
tube to an acceptable value (less than 0.05 inch). A solution which met 
all desired specifications (1.00 £0.03 over . 10 - 100 hertz range) was 
found (see Figure 6) and physically constructed in order to obtain experi- 
mental results to confirm the configuration selected. The configuration 


specifications were: 





Tube Tube Tube 

#1 #2. #3 
Transducer volume (in.3) 0.0 CEO 0.0238 
Tube length (in.) 0.45 29.0 2 V-40 
Tube diameter (in.) 0.0465 0.078 Ore 


The three-tube solutions again confirmed the phase shift prediction; 
the experimental values were within 7% of their predicted values (see 
Figures 8a-d). However, as was the case with the two-tube solution, 
the three-tube experimental configuration produced pressure ratios less 
than those predicted by Bergh's solution for frequencies above 30 hertz 
with the values lower by as much as 30% (see Figures 8a-d). This 
characteristic made the theoretically accentable solution (Figure 8a) of 
limited practical value since the experimental response was consider- 
ably below the prediction. Even though the actual solutions were not 
as predicted, Figure 10 shows trends comparable to those predicted by 
Bergh's solutions for identical configurations (Figure 9). 

From Figures 8b-d, it was observed that changing the lengths of 
tubing raised or lowered the overall curve and a decrease in transducer 
volume (compare Figures 8a, b) tended to increase or raise the overall 


curve in the region of interest and decrease the slope above 80 hertz. 
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A small volume was tested and a possible solution of 1.00 + 0.09 over 

a 10-95 hertz range was finally attained (sige Figure 8d). This configu- 
ration appears acceptable for direct analog data processing in low-speed 
fields of interest. Much of the separated flow phenomenon is involved 
at non-dimensional frequencies or Strouhal numbers, S, less than 0.4 
which would typically correspond to frequencies less than 100 hertz for 


a five-inch diameter circular cylinder at a tunnel speed of 100 feet per 


second. The solution configuration is specified below: 





Tube Tube Tube 

+) + 2 + 3 
Transducer volume (in. 3) 0.0 0.0 OO 74 
Tube length (in.) 0.45 eS) 2120 
imibe diameter Gin.) 0.0465 0.078 O08 


In addition to the tube lengths and transducer volumes, the three 
tube diameters were very important parameters when attaining the theo- 
retical response ecireck The diameter effects and trends were not 
examined, however, due to time and facility limitations and particularly 
due to the limited selection of commercial tube diameters available with 
acceptable inside diameter specifications. 

The purpose for attaining a flat Arexsecines response was to allow 
easy auto- and cross-correlation of a body ina low-speed unsteady 
flow employing a simple, low cost analog method shown by Figure ll. 
The cross-correlation coefficients, which relate the coherence of two 


physical measurements, say at stations (1) and (2), are in general 


defined for a specific value of time lag, oi Ov: 
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<Cyh) Ct ee 


ms (7) : LEP ee <e,*> 2 


However, the elementary analog correlation circuit shown in 
Figure 11 does not have a time lag capability and hence the benefit of 
simplicity is associated with a loss of information. Therefore, the 
information ‘to be obtained must be at zero time lag as specified by 


equation (3): 


R. (o)= Serle) €2 (0) > 
os <ep>2 ge2 5% 
oince the system In Consideration, 1.e., 2 remote pressure trance 
ducer system, can be likened to a typical transfer function system shown 
by Figure la, b, the cross-correlation inputs, ey and Cn: can be repre- 
sented by the transducer outputs, Py and Po: If the disturbance were 
of a sinusoidal nature as it was for the laboratory experiment, then the 


two measuring stations become strictly functions of the frequency, 


Suppose: 


Pi = Po; sin Lap 
P2 = Poxw Sin (wt-@) 


then: 
x, (+) = Po, G lw) Sin (wt- ¢,) 
Br (+) = Por Gv) sin (wt- @- 4,) 
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Substitution into equation (3) will yield a cross-correlation coefficient 
as follows: 


(0) = G,(w) G2 (~) Cos (6+ ¢.- 4 ) 


Riz Gta Ge (2) 


From equation (4) it can be seen that for the case of a pair of information 
signals at a unique circular frequency, &), with the respective transfer 


functions G, (@) and G, (w) equivalent, then: 


i 


Rio (0) = cos @ 


Z 


and a flat frequency response is not necessary. Itis only necessary to 
have identical transfer functions used with both signals. 

However, this analysis eliminates one of the most important types 
of flow fields, the unsteady random signal such as that produced by flow 
separation from a bluff body. In such cases the cross-correlation 
coefficient can no (oneen be simplified to equation (5). Insteadaspecauce 


of its reliance on G as a function of 4 , the output of the cross- 


correlation system becomes: 


~ i G, (w) Gy Cio) g ,, (w) d ty 
Riz (0) ji oo y 
| 1 | J. G,"(w) f, es) deo] | S62) §, 4] 6 * (6) 





where 


p il = power spectral density of signal l 


$ 99 


$ 12 = co-spectral density of signals 1 and 2 


II 


power spectral density of signal 2 
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At the input of the transducer system, the cross-correlation coefficient 


is not a function of G and thus reduces to the form: 


Rie (0) = ___Js_Gntede 
s | f~ G, (wo) deo] fF" Fw) des] ¢ /2, (7) 


In order for (6) to equal (7) for any random signal, the following relation- 


ship is obviously necessary. 


G, (wo) = G, (w) = 1.0 (8) 


Should (8) be the case, then, 


Ry 20) input ~ Ry 210) output 


and the system no longer relies on the gain of the transducer system. 
A complete derivation of the cross-correlation equations may be found 
men ppendix ©, 3 
Because of the flat frequency response, the collected data may be 
inputted directly to the cross-correlation circuit shown in Figure 11 and 
the output monitored by either a D.C. voltmeter or an x-y plotter to 
1 

| . 2.2 DA | 

obtain quantities Ces,» , 6e, » , and ce, » , which allows a simple 


calculation of the cross-correlation coefficient at zero time lag. 
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IV. CONCLUSION AND RECOMMENDATIONS 


Although the breakdown of Bergh's theory seems to be directly 
proportional to the number of tubes in the configuration, it does predict 
fairly accurately the trends that each parameter impresses on a system. 
There is, pereoeee a good possibility of obtaining a flatter frequency 
response ar the one obtained in Figure 8d. Time limitations required . 
the termination of the experiment and the three-tube configuration was in 
no way completely examined. 

It should be noted that the flat frequency response traits considered 
herein applied to a fairly restrictive range of frequencies, namely less 
than 100 hertz. The processing of a broader class of signals, i.e., a 
higher range of frequencies, will probably require the utilization of several 
steps or procedures, including: 

a. Recording of the analog form of the data signals on a magnetic 
tape recorder. 

be Digitizing the analog signals using a common time base or 
elock, and 

Ge Processing the digitized data by means of a digital computer. 

Portions of this procedure have been within the scope of the school's 
experimental facilities within recent times, cf., Jesberg [7] - however, 
new equipment additions or replacements make this former technique no 


longer practical until a repeat is done of the work by Jesberg [7] : 
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It is observed from Figures 7 and 8a-d that the theory under- 
predicts the pressure ratio for both two- and three-tube configurations, 
yet the phase prediction is very close to the correct value. 

The following recommendations are therefore made: 

1) Further experimentation of parameter effects in order to 

obtain a flatter frequency response. 

2.) Possible modification of the computer program to allow 

compensation for the difference patna the actual and 


theoretical response. 


4l 





APPENDIX A 


THEORETICAL RESPONSE OF A 


REMOTELY LOCATED PRESSURE TRANSDUCER SYSTEM 


Bergh's theoretical model applies to any number of transmitting 
tubes in a cascaded tube system as shown by Figure lb; however, many 
assumptions are made in order to reduce the governing equations toa 
workable form. These assumptions are as follows: 

a) the sinusoidal disturbances are very small: 

b) the inner diameter of the tubing is small with 

respect to its length: 

c) laminar flow exists throughout the system: 

d) the unsteady disturbances are perturbations 

imposed on the average values on the parameters 


indicated below: 


1) p= p+ Brad: Pressure 

2) De Po + Ren Density 

3) T = Ts “ Tea Temperature (A.1) 
4) ci US + jen Axial Velocity 

) i vit feet Radial Velocity 


These assumptions allow the governing equations to be reduced to the 


following forms: 
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Navier-Stokes Equations 


ee 12 atu, |) WK )_ 
fo OX +Hpo| SS ty ge [edo (A2) 
- <P =0 (A2a) 
Continuity Equation 
OV V , 
aes ov 1 3 
pe[ Se + Sha 4] 4 WwW, (A3) 


Perfect Gas Law 


RoT 
<= iS ea pP (A4) 


Energy Equation 








Teo . , 
A] Sate Se] = Lop te—.g SPT As) 


Since the equations are now partial differential equations, boundary 
conditions can now be specified: 


a) At r =R, the tube's inside wall, 


1) eer) ea) 
2) 0 
3 T(X,R) = 0 
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b) At or = 0, the tupeseenrer. 


1) v(x,0) = 0 

2) u(X,0) = finite 
3) T(X,0) = finite 
4) p(x,0) = finite 
5) 6 0) ES saboy tits 


From equation (A2a) it is evident that p = p(x), which allows a 


separation of variables technique. 


T= £(x*) h (2) (A6) 


where: 
Lond /2. 
Ve 


pee ex fe = Pranote NumBer 


ve Le R ( 


ia 
R 


Equation (AS) takes on a Bessel equation form of the following type in 


z= Pr 12. 


this type of solution: 


fe ge 
h (2)= Co (@) + C2 Yo (2) + 24cp Flx) (47) 


Apply boundary conditions a3 and b3 


| Pp) 
= aes Serra 
: CyJo(aPr™) Pod Cp (a8) 
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Using (A7) and (A8) in (A6) yields: 


Hx Pilg 
p= <p ee Seal \— To (RX Pr ) (A9) 
ae x Jie (x Py) 


Similar methods using boundary conditions al and b2 can be applied to 


solve equation (1) as given below: 


Lx 
7 Fae ) 7 | dp 
JR (x) a Co Do al x (A10) 


Substituting equations (A9) and (Al10), integrating with respect tor, 
and applying boundary condition a2 will put equation (A7) in the form 


oy ek wet Ta (x Prt) 
F (x) = az o> x ate ¥ S; (% Pr V2) 


_ wd Pp Se) 
2 dxe a as 


Pimce v0) —) Oeimplies Filo .— 0, 


(Tole | ell So Pr'’t) Se Te) dep 
as [ Y Jo (XP 2) i. T, (*) Fee O 





The above equation can be solved for the homogeneous solution 
XK é 
p= A ef + Be $x 


y { 
vy a, | pt To (x PP) I. 
Ms Je (x Pt) 


where y 
wt To *(o) 


9 OF Thm) 
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The coefficients A and B can be determined only after the boundary 
conditions for both ends of the tube have been specified. It should 
also be noted here that the Bessel functions in the equations are complex 
numbers since QK= & (i°/2 y). This is discussed further in Appendix B. 

Two additional assumptions must now be made if the solution is to 
cover an N-tube, N volume configuration: 

Cc) ° the pressure and density in the instrument volumes 

are time dependent only 


d) the pressure expansion in the instrument volume 


is a polytropic process described by 


Py 
P." 


= constant 





The constants A and B are now determined in the following procedure. 
For any tube j, 
at x, = 0 
J 


(A11) 


Ae By 
aio! j j 


Pp= nzetitt 4 Bie Mt 


. 463 £ Tp (*) ie 
Ue = — © gel 
oe Ls) “1h } agers —aer% | i 
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Equations (All and (A12) allow solution for = and B, 
A. = (ee P3-1 ee $5 43 


e az 


$ © Pst3 ~ e~ $545 
B. = Pe @ P3h3 Pa 


etils _ e- Fgls 
The mass leaving the tube can now be calculated: 


2: Gs 2 
m, > 5, Ps, Us 21 Ydr = 





1 Re ds Ty (K) 
yg v é 3 sl. a * 
elt eo proetts ed 
Aw Ta(&; ) Jaye Bye | 


For any tube j + l, 


at xX. al (A13) 


at X44) = L yay 


: = P+ l sei _¢ “+4 L ~ + | 
P sti - A set € + Gens ‘Ss 3 & (A14) 


r ; aoe — 
fo 4) Ds. tl” OF 
Psst Je s 37 
The mass leaving tube j+l is now calculated: 


R- 
LAE fe 0. us, 2dr = 


Psst 





he, Det Tz CX 341) (A 


sth Os, ) 
pare AV; (As +1) 
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Equations (Al3) and (A14) allow solution for Bo and Ba 


~@e,,L+ 
Pst1 - PEC P54 +t 


A... 
at! eo Fsail se am en Pyti bye 


ppcteh! — pss 
e Fatih sa _ e~ Pati lye 


San = 


The mass variation within the instrument volume must equal the differ- 


ence in the mass leaving tube j and entering j+l, namely: 


er at Wi 





Cab (m;,-mz,) e — (Als) 


Substitution of the pressure coefficients into (Al5) and recalling the 


definitions 


eet en wt | evt, o-wt 
cosh wt Ao - and sinh @t = 
a C 
the final pressure ratio prediction is obtained. 
Ps | oor 8 Ws (ayth.) Toles ) 
—— = | Cosy L) + Sa ee 8 Ka we SNH (dL 
Ps (els) TRIO, Tre) 85) 
Z. 
+ DRav1_ Path Tey) Sols) sin (sla) 
WRE 4 Togas) Ta lmy) Sin ($j41 Late) 


: — | 
; cost ($54 L541) — a4 | 
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APPENDIX B 


COMPUTERIZATION OF THE THEORETICAL RESPONSE 


OF A REMOTE PRESSURE TRANSDUCER SYSTEM 


A. Computer Program Development 


In Bergh's theoretical development (Appendix A), the Bessel 
| | | eo 3/2 eee 
function parameter, ©, contains the imaginary term, i , Which can 


, and Bei,. 


easily be represented by the Kelvin functions Ber, Bei, Ber, 


These functions were used in the program in the form of an infinite 
; | -10 
power series and terminated when the terms become smaller than 10 . 


The series expansions were taken from McLachlan [4] and altered into 


the form below: 


Jo(%) = ber yt bel Y (B1) 
y 


+ & We, 
| - y i aan a Y tT os 
(2-4)" (2+ 4-6-8) (2°46 8-10+12)* 


Fake R (&)* = 7 Ve 





a 


O- 
(0 
my 
“ 
iI 


2 G 10 
ber ys D - y Y 


Sennen —_— 
PGE EBB) LEGO EI IO) 
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J, (*)= ber, y + £ bel, y 

-- ye yt y& (B2) 
bey, an aie Lt ae 24.2%, 3 Brice ae 2 

_ Y¥ ye vu y ¢ 
bel, Y > oye ee 24253” Fase °° 





J; 
y= ST 
(B3) 


Bi Computer Program Operation 


The computer program was designed primarily for use on a time- 
sharing system remote terminal where quick output could be obtained 
for a given input. However, it may also be modified for batch processing 
as explained below. The computer program allows selection of either 
al, 2 or 3 tube configuration by the operator and prints out the pressure 
ratio and phase shift for multiples of 5 hertz from 10 to 120 hertz. The 


flowchart is given in Figure (Bl). 
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APPENDIX C 
CROSS-CORRELATION OF A RANDOM UNSTEADY SIGNAL 


The cross-correlation coefficient may be defined in its most 


general form similar to methods described in Bendat and Piersol [8 }. 


<e, (+) C2 (+4+T) » 


R, (+) = 
° Ket) “2 <er) V2 (Cl) 


However, because of the zero time lag limitation of the analog circuit 


shown in Figure 11, this reduces to the form: 


Ce,C2z> 
R 2 =) A a re ee - 
¢¢) Cry 2 eek > le me 


For the system input, the following terms are defined as: 


Kee. = bim 52 J et) e2(4) dt 
_ G , (w) du 
So B Ww)das 


a eS Gs) du 


(C3) 


J] 


<e;? y 
Sen 


}) 
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d 1] = power spectral density of signal 1 


$.. = power spectral density of signal 2 


21, 


After passing through the system, the quantities have a transfer function 


co spectral density of signals 1 and 2 


imposed on their true output, and the following terms are defined: 
nie 4 
E(w) = G (w) e eC, (w) 


e,(w) = A, ror (CA) 


thus: 
<e,ez Vi = (ay G, (w) Go (uw) g., (us) edt 
Ce" > = SS Gy (ws) ¢, (wo) dw 
<e@,> = a Gy Ww) P22 Ww) dw 


A relationship between the output and input cross-correlations may be 


Betained by direct substitution of (C3) and (C5) into (C2), yielding: 


Jo bs) dus 


‘ie G,(w) SG, (a>) € zw) dus 


[So erew) Flu) de [SP et) Fw )do] 2” 


Os: 





There are two possible simplifications available at this point. If the 


input were of a harmonic nature, say ac frequency GW, then the expres- 


sions for spectral density would simplify and yield: 


Sep >=). Fu) Lv-ts,) dao = Fy, (w,) 


Ket aS SF lw) §w-a,) dw = Bar 1) 


(C7) 


Kee. 252 FW) $(w-w,) dw = 2, @) 


similarly: 
< e; >= Pay $1, (w) G,* (as) 5 (us-c,) BGS = $, (t,) GF (uw) 
ie - >= . b , &)E2w) SWw-w,\deo = $,, Ww) SG Ww) 


(C8) 
<6,2,9=56 5,0) GC) Ew) § @s-0) das = $9 E(u SL) 


and the dependence of the zero time lag cross-correlation coefficient 
upon the system's transfer function would disappear due to cancellation. 


However, if the input were of a random nature, then the only 
simplification readily available is to have: 
G, (w) = G 


0 (vs) = US 


in which case the original input cross-correlation is obtained since the 
system transfer function factors out from the integrands, cancels, and 
simplifies to equations: 

<e@;> = ¢e,>= So P (4s) du 


CO2> = <ez>= 55° $0) dw c9) 
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€é,é2) = eHiCcem, = vo 2. (ts) dw 


Thus: 


Riz Co) = R\. (0) (C10) 


Since the latter case also applies to the harmonic signal data 
processing, this is a most desirable situation and is the motivation 


behind seeking a flat frequency of 1.00 + 0.05. 
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